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Abstract 

Fourth-rank tensors [[V 2 ] 2 ] (Voigt's) type, that embody the elastic properties of crystalline 
anisotropic substances, were constructed for all 2D crystal systems. Using them we obtained 
explicit expressions for inverse of Young's modulus E(n), inverse of shear modulus G(m, n) and 
Poisson's ratio z/(m, n), which depend on components of the elastic compliances tensor S, on di- 
rection cosines of vectors n of uniaxial load and the vector m of lateral strain with crystalline 
symmetry axes. All 2D crystal systems are considered. Such representation yields decomposition 
of the above elastic characteristics to isotropic and anisotropic parts. Expressions for Poisson's co- 
efficient are well suited for studying the property of auxeticity and anisotropy 2D crystals. Phase 
velocities are calculated for all 2D crystal classes. 
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I. INTRODUCTION 



In our previous paper we derived explicit expressions for inverse of Young's modulus 
E(n), inverse of shear modulus G(m, n) and Poisson's ratio z/(m, n) for 3D crystal systems 
of high and middle symmetry Our derivation yields expressions which do not depend 
on the choice of Cartesian coordinate system. These expressions depend on components of 
the compliances tensor S, and on direction cosines of vectors n of uniaxial load and the 
vector m of lateral strain with crystalline symmetry axes. To obtain such formulas, we used 
fourth-rank tensorial bases introduced by Walpole j2j. He considered such tensorial bases for 
all thirty-two 3D crystal classes. The fourth-rank tensors of Voigt's type for 2D crystalline 
solids do not seem to have been submitted previously, therefore, we construct and use them 
to calculate the above mentioned characteristics of elastic media belonging to all four 2D 
crystal systems 

Auxetic elastic materials exhibits unusual property of becoming wider when stretched and 
thinner when compressed, i.e. they exhibit a negative Poisson ratio. Love |4| presents a single 
example of cubic single crystal pyrite as having a Poisson's ratio of -0.14. Lakes described 
the synthesis of an actual auxetic material and proposed a simple mechanism underlying the 
negative Poissons ratio Alderson and Evans fabricated microporous polyethylene having 
a negative Poissons ratio (|. Data were compiled from the literature by Baughman et al. 
and were analyzed to show negative Poissons ratio to occur for stress in an oblique direction 
upon single crystals of cubic metals pj. 

The unusual properties of auxetics imparts many beneficial effects on the material's 
macroscopic properties that make auxetics superior to conventional materials in many com- 
mercial applications. To explain microscopic mechanisms underlying the property of aux- 
eticity, various models are proposed. Two-dimensional microscopic crystalline models make 



an interesting class of auxetics (cf. e.g. 8(- Each model of such kind belongs to one of 
crystal systems, therefore its auxetic properties are characterized by Poisson's ratio suitable 
for this system. 

For completeness, in the Appendix we consider the Christoffel tensor for all 2D crystal 
Laue groups in the Appendix to this paper. 
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II. BASES OF FOURTH-RANK TENSORS FOR 2D CRYSTALLINE MEDIA 



Repetition of a latin suffix in a product of tensors or in a single tensor implies summation 
with respect to that suffix over values 1,2. If u, v are two mutually perpendicular vectors, 
then 

UiUi = 1, UiVi = 0, UiUj + ViVj = Sij, (I) 

the last of these relations can be written in a dyadic form (u <g> v + v <g> u) = 1%. 
The stiffness tensor C and the compliances tensor S are mutually inverse, i.e. 

SC = CS = h. (2) 

The product AB of two fourth-rank tensors A and B has components (AB)jki = Aij rs B rs ki. 
The fourth rank identity tensor J 4 has components 

(h) ijkl = (5 ik 5ji + SaSjk) /2. 
A. Hexagonal (isotropic) system 

Consider crystal belonging to the hexagonal (isotropic) system. Two tensors J and K 

JijM = SijSki/2, Kijki = (5ik5ji + 5u5jk — SijSfd) /2, 

with the multiplication table J 2 = J, K 2 = K, JK = KJ = bring about the decompo- 
sition J + K = J 4 . The general fourth-rank tensor of complete Voigt's symmetry (i.e. of 
[[V 2 ] 2 ] type) is constructed as 

A = ajJ + axK, (3) 

where 

a j = Aujj/2, a K = {A m - aj) /2. (4) 
Applying Eqs. © and (jlj to C and S tensors, we obtain 

cj = Cxi + Cx2, cr — Cxi ~ Cx2, 
sj = Sxx + S12, Sjc = Six — S12, (5) 

where C n = Cxxu, C 12 = C llw . Similarly, S n = S ni i, S 12 = S U 22- 
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An isotropic 2D medium is mechanically stable if 

sj = S n + Si 2 > 0, s K = S u - S 12 > 0, (6) 

i.e. Sn > \Si 2 \. From inequalities (jUJ) it is easily shown that Sn > 0. Eigenvalues cj, ck 
and Cn, C12 obey analogous inequalities. 

Using Eqs. (J2J) and multiplication table for J and K, we express components of S by 
components of C 



C u n C 



>^2 ' 12 r>% ' ^ ' 

°11 — °12 °11 — °12 

B. Quadratic system 

Two crystallographic directions of the quadratic system, a and b, are mutually perpen- 
dicular, hence they obey the relations (JTJ. Now the tensors K and 1± decompose 

K = L + M, I A = J + L + M. 

The tensors L = (u Cg> u — v ® v) ® (u ® u — v ® v) /2 and M = (u ® v + v £g> u) £g> 
(u ® v + v ® u) /2 are idempotent and have components 

£yfcj = (uiUj - ViVj) (u k ui - v k vi) /2, (8) 
Mjfez = {uiVj + UjUj) (MfeUj + v k ui) /2. (9) 

The tensors J, L, and M fulfill the multiplication table 

L 2 = L, M 2 = M, JL = LJ = 0, JM = MJ = LM = ML = 0. (10) 

The general fourth-rank tensor of the [[V 2 ] 2 ] type is composed as 

A = ajj + a L L + a M M, (11) 

where 

a L = (uiUj - ViVj) A ijk i {u k ui - v k vi) /2, 

a M = (uiVj + ViUj) Aij k i (u k vi + v k ui) jl. (12) 

Using Eqs. (fTT|) and (|12|) for C and S 1 tensors, we obtain in the reference axes coincident 
with the crystal reference axes (CRA for short) 

Cl = (Cll — C12) ) C M = 2Ci212, S£ = (Su — ^12) , % = 25*1212. (13) 
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A quadratic medium is mechanically stable if c a > or s a > (a = J,L,M). In terms of 
S T j J,J = 1,2,6 

S n > 0, S m > 0, (14) 

and conditions (jH} are fulfilled, i.e. Sn > \Si 2 \. 

The components of S are related to components of C. Eqs. (0) are fulfilled and, since 
C m = C1212 and S 66 = 4S , m 2 Q, 

5 66 = 1/C 66 . (15) 

When ai = (Im, the general symmetric fourth rank-tensor A reduces to the corresponding 
tensor for the isotropic 2D system. For A = C, one obtains the familiar relation Cqq = 
(Cn - C12) /2 and S m = 2 (S u - S 12 ) for A = S. 



C. Rectangular system 

As the crystallographic directions of all classes of the rectangular system are mutually 
perpendicular, the conditions (JTJ) are satisfied by the unit vectors u, v that define these 
directions. The tensor M is retained intact. The four new tensors are introduced 

E u = (u ® u) ® (u ® u), E u = (u ® u) ® (v ® v), 
E 2 i = (v ® v) ® (u® u), E 22 = (v ® v) ® (v ® v). (16) 

The tensors M, £?p ? , (p, g = 1, 2) form the complete set whose multiplication table is The 
general fourth-rank tensor of the [[V 2 ] 2 ] type is constructed by the linear combination 

A = a u E n + a 22 E 2 2 + a s E s + a M M, (17) 
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where E s = (E12 + E21) and 

a n = UiUjAijkiUkUu a 22 = ViVjA ijk iv k Vi, 
a l2 = UiUjA ijkl v k vi = a 2 i = a s , 
a M = (ujVj + ViUj) A ijkl (u k v t + v k ui) /2. (18) 

In CRA, by combining Eqs. (|T7|) and (fTHj) for C and 5* tensors, we can evaluate c Q/ g (a, (3 = 
1,2, M) as follows 

C\l = Cxi, C22 = C22, C12 = C21 = C12, Cm = 2C*i212 = 2C*66, 

«ii = 5*11, S22 = S22, S12 = s 2 i = 5*12, s M = 25*1212 = 5*66/2. (19) 

It follows from Eqs. dHJ) and (HHJ) that 

J = (E n + E22 + E s ) /2, K = (E n + £ 22 - E s ) /2 + M, 

L = (E n + £ 22 - E„) /2, J 4 = £11 + £22 + E M . (20) 

In the CRA, relations connecting compliances and stiffnesses of the rectangular system 
read 

_ C22 „ _ Cn 

W1^22 — ^12 — W 
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C12 

5i2 = JyTi *S*66 = Cq 6 . (21) 

<^ll<^22 — ^12 

A rectangular elastic medium is mechanically stable if inequalities (|14|) are fulfilled and 

522 > 0, 5n522 > 5^2- (22) 

D. Oblique crystal system 

Two unit vectors a and b make an arbitrary angle (7^ tt/2) with each other in the 
plane that they define. It is convenient to introduce two further unit vectors u, v which 
are mutually perpendicular as independent combinations of a and b, e.g. a = u, v = 
[b-(ab)a]/[l-(ab) 2 ] 1/2 Q. 
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The tensors E aa (a = 1,2) and M = E 33 are retained intact. Four new tensors are 
introduced 

E13 = (u ® u) <S> (u <g> v + v <S> u) /a/2, 
£31 = (u <g> v + v <g> u) ® (u <g> u)/\/2, 

^23 = (V <g> V) <g> (U <g> V + V <g> U) /\/2, 

£32 = (u <g> v + v <g> u) <g> (v <g> v) /\/2. (23) 

The multiplication table 
can be now compiled. 

The pairs tensors E 12 and E 2 \, E n and E 3i , E 2 3 and E 32 have to be made symmetric 
if the general symmetric tensor is required, hence we introduce three symmetric tensors, 
namely 

E§ = (E aj3 + E Pa ) (a ^ (3, a, (3 = 1, 2, 3) (24) 

Using tensors E aa and E^l (a 7^ (3), (a, (3 — 1, 2, 3), the general fourth-rank tensor of [[V 2 ] 2 ] 
can be constructed for the crystal 2 system 



where a a p = a^ a for a 7^ (5 and 



Ajki - ^2 aa P E ijkii 



(25) 



(26) 
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When ai 3 = a 23 = 0, Eq. (JHJ) reduces to Eq. (fTTjl. 

The formulas ()25j) and bring about decomposition 



J 



+ -^22 + E 



/2, K 



2:2 



/2 + E 



33 • 



-En + E 22 — £?i2 



/2, 7 4 = E n + B 



22 T £/33 



(27) 



Invoking the definitions of the tensors E a p and Eq. (|2fi|) . we can calculate the elements 
of the matrix representing the tensors C and S in the CRA 



Cll — Cii, C 2 2 — C22, C33 — 2C 66 , C12 — c 2 i — C\2, 
Cl3 = C31 = v2Cni2 = V2Ci6, C 23 = C32 = V / 2C 2 212 = V2C26, 



(28) 



Sll — Six, s 22 — S22, S33 — Sqq/2, Sx2 — S21 — 5i2, 

sx3 = s 3 i = v / 2S' 1 ii2 = Sxe/V2, s 23 = s 32 = V2S2212 = S2&/V2. (29) 

Using Eq. (|2"|). we obtain the set of six linear equations for s a/ g. Solving them we obtain 
expressions relating the elements of the matrices representing tensors S and C in the CRA 



Sxi 



< ~ / 22 ( -'66 — ( - / 26 n ^12^66 — ( - / 16 L/ 26 
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D 



_ 022^16 — ^12^26 c _ ^11^26 — ^12^16 

where D = (C11C22C33 + 2C12C26C16 — C22C16 — C11C26 ~~ C\ 2 C^. 
The stability conditions for oblique crystalline media read 



•Sii > 0, S22 > 0, Sqq > 0, S'iiS'22 > S" 



12; 



S'll'S'26 — 2S'i25'265'l6 + 5 , 22>S'l6 < ^ (*S'll*S'22 — S^) ■ 



(30) 



(31) 



III. ELASTIC MATERIAL CONSTANTS FOR 2D CRYSTALLINE MEDIA 



Consider any two specified orthogonal unit vectors n and m and four elastic material 
constants of a considered body, namely the bulk modulus k, Young's and shear modules 
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E(n), and G(m, n) respectively, and Poisson's ratio z/(m, n). They are defined as 12L Il3| 

k- 1 = I 2 ■ S ■ I 2 , (32) 
(n <g> n) • S ■ (n <g> n) = nin,jS ijM n k n u (33) 



£(n) 

(m <g> n) • S 1 • (m <g> n) = minjSij k im k ni, (34) 



4G(m, n) 

i>(m, n) . . 
; = (m ® m) • 5 • (n <g> n) = mimjb ijk in k ni. (35) 

Since k, i? _1 (n) and [4G(m, n)] -1 are defined by quadratic forms of positive definite 
tensor, they are positive, i.e. k > 0, and £' _1 (n) > 0, [4G(m, n)] _1 > for all mutually 
perpendicular vectors n and m. If u(m, n) is negative an elastic medium is said to be an 
auxetic. 



A. Isotropic 2D media 



Using definitions (|35| ) -()35| ) . we obtain 

E 1 = sn = S n = —2 — ^2"> p = ~ s i2 = -S12 = — ^2~' ( ' J)(i) 



G- 1 = 2s K = 2 «£i = 2 SJ = 2 (5 n + S 12 ) = ^ ■ (37) 

All these elastic material constants are isotropic, and when conditions (JHJ) are fulfilled, k, E 
and G are positive. For all remaining crystal classes the bulk modulus is equal to K is . 



B. Quadratic 2D media 

In the case of 2D quadratic media E(n), G(m,n) and z/(m, n) are anisotropic, as they 
depend on the directional cosines (nu), (nv) and (mu), (mv) 

E-\n) = sj/2 + s L [(nu) 2 - (nv) 2 ] 2 /2 + 2 SM (nu) 2 (nv) 2 , 
[4G(m, n)] _1 = Sl [(mu)(nu) — (mv)(nv)] 2 /2 
+s M [(mu)(nv) + (mv)(nu)] 2 /2, 
-z/(m, n)/E(n) =sj/2 + s L [(mu) 2 - (mv) 2 ] [(nu) 2 - (nv) 2 ] /2 

+2s if (mu) (mv) (nu) (nv) . (38) 

For mechanically stable media, both E(n) and G(n, m) are positive for all mutually 
perpendicular vectors n and m. 
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C. Rectangular 2D media 



In the case of rectangular 2D crystalline media, expressions for elastic material constants 
are more complex than for isotropic and quadratic media 

E-\n) = su(nu) 4 + s 22 (nv) 4 + 2 (s 12 + s M ) (nu) 2 (nv) 2 , 
-z/(m, n)/E(n) = Sn(mu) 2 (nu) 2 + s 22 (mv) 2 (nv) 2 
+si 2 [(mu) 2 (nv) 2 + (mv) 2 (nu) 2 ] + 2%(mu)(mv)(nu)(nv), (39) 
[4G(m, n)]- 1 = Sll (mu) 2 (nu) 2 + s 22 (mv) 2 (nv) 2 
+2si2 (mu) (nu) (mv) (nu) + sm [(mu)(nv) + (mv)(nu)] 2 /2. 

D. Oblique 2D media 

In the case of oblique 2D crystalline media, we obtain 

E-\n) = sn(nu) 4 + s 22 (nv) 4 + 2 (s 33 + s 12 ) (nu) 2 (nv) 2 
+2 v / 2si 3 (nu) 3 (nv) + 2 V / 2s 23 (nu)(nv) 3 , 
— z/(m, n)/E(n) = Sn(mu) 2 (nu) 2 + s 22 (mv) 2 (nv) 2 
+2s 33 (mu)(mv)(nu)(nv) + s 12 [(mu) 2 (nv) 2 + (mv) 2 (nu) 2 ] 
+^13 [(mu) 2 (nu)(nv) + (nu) 2 (mu)(mv)] 
+^^23 [(mv) 2 (nu)(nv) + (nv) 2 (mu)(mv)] , (40) 

[^n)]- 1 = Sll (mu) 2 (nu) 2 + %(mv) 2 (nv) 2 
+s 33 [(mu)(nv) + (mv)(nu)] 2 /2 + 2si 2 (mu)(nu)(mv)(nv) 
+v / 2 [si 3 (mu)(nu) + s 23 (mv)(nv)] [(mu)(nv) + (mv)(nu)] . (41) 

It should be noticed that vectors of stretch n and strain m are two mutually perpendic- 
ular vectors of the unit length. In the reference frame related to the symmetry axes their 
components fulfill three conditions 

(nu) 2 + (nv) 2 = 1, (mu) 2 + (mv) 2 = 1, (nu)(mu) + (nv)(mv) = 0. (42) 
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This means that only one component of n and m is independent. That implies, via 
that mi = —ri2 = — sin tp, and m 2 — n\ — cos<£>. The explicit expressions for elastic material 
constants of 2D crystalline media and their consequences are studied in a companion paper 

0. 



Appendix: the Christoffel tensor 

The Christoffel (sound propagation) tensor T(n) is defined as jl^ 

r(n) = n • C ■ n. 

The Christoffel equation allows one to find the phase velocities c(n) and polarization 
vectors U of a crystalline medium with the density p 

(T rs - pc 2 5 rs ) U s = 0. (43) 

With the help of Eq. (J25J) for oblique system, we obtain 

r o (n) =7 mi (n)(u®u) + 7„„(n)(v <g> v) + 7 u „(n)(u (g) v + v <g> u), 



where 
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(o), 



n) = c n (nu) 2 + c 33 (nv) 2 /2 + v / 2ci 3 (nu)(nv) 



7 ^(n) = c 22 (nv) 2 + c 33 (nu) 2 /2 + v / 2c 23 (nu)(nv), 
7^ } (n) = 7^(n) = (c 12 + c 33 /2) (nu)(nv) 
+Ci 3 (nu) 2 /v / 2 + c 23 (nv)7V2. 

For rectangular, cubic and isotropic systems, Ci 3 = c 23 = 0. For quadratic system 
Cn = c 22 , and for isotropic media cl = cm- 
In CRA one obtains 

7|° } (n) = C u nl + C m nl + C x& n x n 2 , 
72? ( n ) = C 2 2n\ + C m n\ + C 26 ?2in 2 , 
7i2 ) ( n ) = 72? ( n ) = (Cia + C m ) n x n 2 + C m n\ + C 26 n 2 . 

Taking into account that for rectangular and quadratic systems C\§ = C 2 6 = 0, and addition- 
ally that for quadratic system Cn = C 22 , one obtains Christoffel's matrix for these systems. 
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If one takes Cqq = (Cn — C12) /2, then one obtains Christoffel's matrix for isotropic media. 
On the other hand, using Eq. (j3J) one obtains 



r«(n) 



cjn\ + c K cjnin 2 
cjn x n 2 cjn\ + c K 



with cj and c^- given by Eqs. (J3J). Christoffel's tensors of isotropic 2D media obtained on 
the above two ways are identical. 

The propagation vector n = and the vector q = XJ t perpendicular to n are the 
polarization vectors. The corresponding phase velocities are 



q = y/(Cj + C K ) ]2p = y/Cn/p, 
ct = y/C K /2p = v/(C n - C 12 ) /2p. 

The Christoffel tensor and phase velocities for quadratic system are discussed in a separate 
paper 
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